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1.0  PURPOSE 


The  purpose  of  this  project  shall  be  the  construction  of  a 
feasibility  model  demonstrating  the  principles  of  optimal  discontin¬ 
uous  control  for  servo  motors.  The  particular  servomechanism  to  be 
constructed  will  be  optimally  designed  with  respect  to  inputs  of  the 
type,  step  plus  ramp  disturbance,  and  will  be  self  adaptive  to  torque 
variations  in  the  load.  There  appear  to  be  two  possible  design 
approaches  to  the  synthesis  of  the  optimal  switching  function.  One 
of  these  would  be  to  construct  an  analog  of  the  switching  function 
using  conventional  analog  computer  techn.lques,  such  as  logarithmic 
amplifiers,  multipliers,  dividers  and  other  nonlinear  elements  to 
produce,  in  real  time,  subject  to  the  state  variations  of  the  error  and 
error  rate  in  the  phase-plane,  an  optimal  discontinuous  forcing 
function. 

A  second  synthesis  technique  would  be  to  construct  the  switching 
function  faster  than  real  time  through  the  use  of  passive  analog 
elements.  This  prediction  technique  would  run  the  linear  plant  to 
the  origin,  in  the  phase-plane,  many  times  faster  than  real  time. 

Noting  the  error  and  switching  the  forcing  function  will  match  that 
predicted  condition.  This  second  technique  has  been  adopted  by  the 
project  managers  of  this  contract  for  the  following  reasons: 

1.  It  appears  to  be  economically  much  more  feasible  than  the 
fabrication  of  many  expensive  nonlinear  computing  elements. 

2.  It  is  expected  that  the  accuracy  of  the  passive  elements 
involved  will  be  higher  than  those  of  the  equivalent  analog 
synthesis  techniques. 

The  technique  appears  to  be  new  and  therefore  represents 
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3. 


a  step  forward.  In  the  state  of  the  control  system's  art 
wif^h  respect  to  relay  actuated  control  systems. 

4.  It  is  not  felt  that  this  technique  is  in  any  real  sense 
a  deviation  from  the  statement  of  work  of  the  contract 
in  the  sense  that  the  optimal  switching  curve  is  indeed 
fabricated,  faster  than  real  time,  and  sampling  techniques 
provide  us  with  the.  answer  as  to  whether  or  not  the  linear 
plant  is  following  the  optimal  switching  function. 

Wliile  on  the  face  of  it  there  appears  to  be  a  loss  of  generality 
in  constructing  a  faster  than  real  time  analog  for  the  second  order 
servomechanism,  it  is  also  true  that  in  any  generalized  n-th  order 
control  system's  problem,  the  use  of  logarithmic  amplifiers  in  ana¬ 
log  synthesis  is  mandatory.  Once  the  system  has  been  reduced  to  the 
n-1  order,  in  phase-space,  the  final  .solution  could  be  mechanized  in 
the  same  manner  that  it  is  being  produced  in  this  particular  project. 

A  completely  general  approach  would  be  to  use  the  logarithmic 
amplifier  synthesis  technique  for  the  entire  control  system's  prob¬ 
lem.  However,  due  to  the  unavailability  of  sufficiently  accurate 
log  amplifiers,  it  is  felt  that  a  feasibility  determination  can  be  made 
on  the  basis  of  this  technique.  Further  studies  will  improve  the 
techniques  required  to  mechanize  the  completely  general  n-th  order 
problem. 
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2 . 0  ABSTRACT 


Using  results  obtained  by  Dr.  R.W..  Bass,  who  introduced  the 
concept  of  relay  feedback  into  the  field  of  control  theory,  it  has 
been  possible  to  mathematically  describe  a  phase  plane  portrait  of 
the  switching  characteristic  of  an  extremal  torque,  time  optimal, 
positional-servomechanism  which  is  self  adaptive  to  torque  load 
variations . 

This  problem  has  been  of  historical  interest  to  workers  in  the 
field,  and,  at  the  present  time,  of  accentiiated  practical  interest 
in  the  field  of  guidance  and  control.  Because  of  this  new  emphasis 
on  extremal  and  optimal  control  the  current  feasibility  effort  was 
instituted.  The  simplifying  assumptions  and  mathematical  idealiza¬ 
tions  have  been  dropped  and  an  engineering  attack  has  been  mounted 
on  the  practical,  the  realizable,  the  economic,  and  the  uncertain 
factors. 

The  approach  has  been  twofold;  theoretical  work  which  removes 
the  original  simplifications,  introduces  uncertain  parameters  and 
considers  the  inherent  difficulties  of  extremal  systems  as  inmedlste 
problem  objectives;  practical  implementation,  as  the  second  major 
effort,  has  been  devoted  to  consideration  of  the  factors  influencing 
the  choice  of  design,  generation  of  system  block  diagrams  (both 
functional  and  operational)  and  the  construction  of  feasibility 
hardware  to  justify  our  choice  of  system  elements. 

Both  of  these  efforts  have  produced  problem  areas  and  compromise 
structures.  In  many  cases  the  choice  of  solution  is  either  arbitrary 
or  must  be  based  on  inadeqvtate  engineering  data  or  theoretical  re¬ 
search  which  has  not  yet  been  accomplished.  In  view  of  these  prob- 
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lem  areas  It  la  felt  that  several  meetings  with  Mr.  Carl  Allen, 
the  Contracting  Officer's  Designated  Representative,  would  be 
advisable  to  resolve  these  Issues. 
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3.0  CONFERENCES  AND  PUBLICATIONS 


3. 1  Conferences 

Dr.  R.W.  Bass,  principal  Investigator  and  W.S.  Hare, 
project  manager,  attended  the  Joint  Automatic  Control  Conference, 
1961,  with  Carl  Allen  of  the  Signal  Corps  and  had  several  important 
informal  discussions  with  leaders  in  the  field  of  optimal  and 
adaptive  control  on  subjects  relating  closely  to  this  project. 

Dr.  Bass  and  Mr.  Hare  met  on  16  October  1961  with 
Samuel  Adler,  Carl  Allen  and  Peter  Devreoctes  of  the  Signal  Corps 
to  discuss  the  status  of  the  project  and  plan  the  future  work. 

3.2  Publications 

No  publications  have  been  approached  to  date,  however 
there  appears  to  be  a  class  of  unusual  problems  represented  in  the 
multipll-connected  phase>space  difficulty  which  would  Justify 
publication  in  advanced  journals. 
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4.0  FACTUAL  DATA 


4. 1  Theoretical  Work 

The  theoretical  work  has  been  devoted  to  mathematical 
effort  to  obtain  a  rigorous  switching  characteristic  in  closed 
form  which  at  the  same  time  includes  consideration  of  the  fact 
that  any  desired  final  position  of  a  rotating  servomechanism  may 
be  taken  modulo  217.  This  means  that  we  have  a  multipli  connected 
phase  space  and  this  is  a  new  consideration  in  this  field. 

In  addition,  computer  simulation  studies  have  been  instituted, 
prlmarly  as  a  means  of  establishing  practical  tolerances  on  the 
system  hardware  to  attain  an  acceptable  limit  cycle,  stability 
margins  and  the  effect  of  stochastic  variations  on  the  input 
signal  and  other  system  parameters. 

In  our  effort  to  derive  a  switching  characteristic  applicable 
to  practical  transducers  in  a  medium  power  servomechanism,  the 
following  calculations  have  been  made. 

Consider  the  following  configuration  for  motor  and  load.  In 
which  the  following  idealizations  have  been  assumed: 

a.  Motor  transfer  function  is  linear  and  second  order. 

b.  Gear  train  is  free  of  backlash. 

c.  Load  parameters  are  constant. 

d.  External  torque  variations  never  exceed  the  stall  torque 
of  the  motor. 

See  Figure  1. 
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Figure  1 

Pictorial  of  Motor  end  Loed  Configuration 
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If  the  servo  torque- speed  curves  have  the  following  form: 


Figure  2 


end  E  >  E2  > 
and  since 


then  let 
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Hence  the  developed  torque 

(1) 


fm  -  Ke  hn  *  »' 


The  load  torque  at  the  motor  shaft 


(2) 


+  T  I  Q  *  + 

1? 


-  r. 


Equating  (1)  and  (2)  we  have  the  system  differential  equations 


Let 


*  J, 


T 


d 


\n  I 

0'  -  n0 


f 


— —  =  F  sgntf"  {the  extremal  constraint} 

Hence  (3)  becomes 

(4)  sJ  0  +  f  0  “  F  sgno^  -I-  Tj, 

If  the  solution  of  this  equation  is  to  follow  a  given  Input, 
0£(t)  so  chosen  that 

F+Td 

(4a)  *  a  +  bt,  /b/  <  — j —  (runaway  velocity) 

within  any  finite  Interval,  then  it  is  convenient  to  define  a  new 
variable,  e(t),  the  actuating  error 


(5)  e(t)  -  (ej^(t)  -  0(t)) 

clearly  |e(t)J-^0  as  0(t)  -^0£(t)« 

After  substitution  of  (5)  and  (4a),  (4)  becomes 

(6)  J  e  +  f  e  -  (fb  -  T^)  -  p  sgn  o' 

We  then  desire  that  choice  of  ^  which  will  bring  e(t)  and  its 
derivative  to  the  origin  from  any  initial  position  (e^,  e^). 

-  8  - 


<  ^  ^ 


Let  us  solve  (6) 

Set  e(0)  * 

^(0)  *  e 

o 

Now  ^ 

e(t)  -  A  +  B£  ^  j^(fb  -  Tjj)  -  F  sgntrjt 

e(t)  “  -  I  B6  ^  Ijcfb  -  T^)  »  F  sgnO"] 

then 

-  A  +  B  _ 

*o  '  •  J  ®  +  T  '  ’'d>  '  ''  ‘^] 

hence 

B  -  ^  |^(fb  -  Tj)  -  F  sgniy^  -  ^ 

A  -  Cq  +  4  ^  [(fb  -  Tj)  -  F  sgn 

Finally 

(7)  e(t)  ■  +  4  ®o  “  ^  ^  sgno-J 

■) 

+  jp  [<^-Td)  ’  ^  +  T^^-Td)  ■  *■ 

and  ^  f 

(8)  e(t)  *  j  [(fb-T^)  -  F  sgn^-J-  -f  ]( fb-T^)  -  F  sgii*-  -  fej6"^ 

If  we  eliminate  the  parameter,  t,  the  resultant  equation  de¬ 
fines  two  families  of  phase  plane  trajectories,  where  the  phase 
plane  is  defined  to  be  e(t),  e(t)  plane.  These  two  families 
differ  in  the  assumed  algebraic  sign  <y  ,  See  the  following  sketches 
wherein  it  is  assumed  that  b,  the  input  rate,  and  T^,  the  disturbing 
torque,  are  held  constant  throughout  the  trajectory 
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-pUn«  Saai-trajcctorlct 


Note  that  only  that  part  of  the  trajectory  which  is  heavily 
shaded  in  the  sketch,  will  bring  the  system  to  the  origin  with  no 
torque  reversal,.  This  then  is  the  desired  final  trajectory.  Let 
us  solve  equations  (7)  and  (8)  for  this  trajectory  by  eliminating 
t,  and  setting  the  endpoint  equal  to  the  origin,  e(t)  ■  e(t)  ■  0. 
Hence,  (8)  becomes 


Clearing  (7)  and  (8)  we  have 

<10) 


-j  /Ufb-Tjj)  -  F  sgn<>‘J  -  fe, 
F  3gn<rJ  log^- - . -"F  's^tfJ 


Equation  (JD)  is  a  double  valued  implicit  function  of  (e^, 

See  illustration  on  page  10a. 

Note  that  In  the  I  and  HI  quadrant  solutions,  it  is  Implied 
by  equation  (9a)  that  we  reach  the  origin  in  negative  time,  l.e., 
an  unrealizable  solution. 

In  other  words,  we  require  that  equation  (10)  holds  in  reality 
only  if 

(11)  sgn«"  -  sgn  e 


Hence 

(12)  0  -  e^  +  ~  ^ 


F  sgn  ejlog  |^1  +^(T^_fbT 


f  sgn  h 


3J 
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Figure  4. 


Second  Order  Switching  Boundary 


.  / 


10a 
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becomes  the  optimal  trajectory,  and  Indeed  it  will  be  shown,  Is  the 
switching  curve. 

If  in  the  derivation  of  (12)  we  let  the  end  point  be  given 

e(t)  =  e^ 
eft)  -  0 

then 

(12a)  +  F  sgn  eJlog|l  + 

which  implies  that  the  system  comes  to  rest  with  an  overshoot  of  e^^. 

Note  that  the  sign  of  e^  is  also  the  sign  of  the  applied  torque  which 
will  cause  the  system  to  evolve  to  a  point  of  interception  with  the 
ideal  trajectory  (12). 

Hence,  if  we  set  ^  and  the  applied  torque  eqt^l  to,  F  sgn  0^ 

then  we  generate  the  complete  closed  loop  ideal  switching  characteristic 
for  arbitrary  initial  conditions. 

(13)  cr-  ^  ^  ^  ■*' 

For  step  response  only,  with  no  torque  adaptation  we  have 

(U)  «■-  *0  +  4  ®o  -  ^  (“*"  ®o) 

See  Illustration  on  following  page,  11a. 

Equations  (13)  and  (14)  are  completely  rigorous  with  respect  to 
the  assvimptions  listed  at  the  beginning  of  this  section;  however, 
we  note  the  sensors  of  a  rotating  servomechanism  are  really 
measuring  the  variables  of  a  congruence  class,  namely 

0^  “  9  +  ZTn  (n  ■  any  Integer) 
is  a  solution  of  the  equations. 

Representing  9  s  0  mod  lH  is  equivalent  to  mapping  the  phase  plane 
on  a  cylinder.  See  pages  lib,  11c,  lid. 
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Figure  5. 

Phase-plane  Trajectories 
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Figure  6. 

Mep' of  Anblguous  Phase-plane. 
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Figure  7. 

Truncated  Phase-plane  Plot 
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Figure  8. 

Cylindrical  Map  of  Phase-plane  (ambiguity  removed). 
(The  'shaded  regions  represent  o^O) 
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The  construction  of  a  computer  analog  of  this  switching  criterion 
is  non- trivial.  Several  structures  present  themselves,  however. 

Perhaps  the  most  generally  useful  method  is  the  analog  s)nithesis. 

This  difficulty  with  the  multipli-connected  phase  surface  la 
avoided  in  part  by  the  selection  of  components. 

For  example  we  note  that  the  motor  and  load  are  in  a  closed 
system  in  which  it  is  assumed 

|Tj  -  fb|  :£F. 

Clearly  then  /e/  ^  is  the  worst  case,  condition.  If  the  e  axis 
is  expanded  through  choice  of  gear  ration  then  the  semi- infinite  strip 

of  Figure  9  (page  12a)  can  be  represented  in  a  finite  box  with  as 

* 

few  as  two  discontinuities  in  traversing  the  range  of  e  or  e.  See 
Figure  10  on  page  l2b. 

Another  source  of  difficulty  arises  from  the  fact  that  no 
synthesis  will  be  perfectly  realized.  Analogous  to  this  are  the 
well  known  oversimplifications  of  the  signum  operator,  the  inherent 
nonlinearities  of  the  plant,  and  the  Inaccuracies  of  the  sensors. 

These  realities  lead  one  to  believe,  that  the  system  will  have 
several  torque  transitions  in  its  evolution  to  the  origin.  Usually 
this  multiple  switching  Is  a  limit  cycle  or  a  chattering  regime  or  both. 

We  have,  in  our  analog  synthesis  program,  observed  these  phenomena 
in  an  attempt  to  put  practical  limitations  on  the  acceptable  tolerances 
of  the  servo  components.  Because  our  multiplier  cabinets  have  been 
inoperative  since  the  outset  of  this  program,  we  have  not  as  yet  closed 
the  loop  on  the  analog  synthesis  of  equation  (13).  We  have,  however, 
approximated  equation  (14)  with  straight  line  segments. 
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Figure  9. 

Certheslan  Plot  of  Phaae-plane  Showing  High 
Gear  Ratio. 
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Figure  10 


Cylindrical  Plot  of  Preceding  Page  Sho«d.ng 
'High  Gear  Ratio. 
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4 » 1 , 1  List  of  Symbols 


T  “  torqxie 

*  disturbance  torque  at  motor  shaft 

*  disturbance  torque  at  load 

Tjjj  ®  developed  torque. 

Tj^  =  load  torque. 

Jj^  =  inertia  of  load 

J_  *  inertia  of  armature 
m 

J  =  equivalent  inertia  at  load 
f  =  load  friction 

f  =  motor  friction 
m 

f  “  equivalent  friction  at  load 

n  =•  gear  ratio 

Kg  ■  motor  torque  constant 

*  motor  damping  constant 
0'  -  shaft  displacement 

0  -  load  displacement 

terminal  voltage 

<3^  -  switching  function  (or  control  function) 

F  -  ^e^in 
n 

*  input  command 
e(t)  “  error 

a  «  initial  displacement  of  input  comnand 

b  “rate  of  change,  of  0£(t) 

=•  initial  error 
o 

“  initial  error  rate 

A,  B  undetermined  coefficients  of  differential  equations 
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£  >  base  of  natural  logarithms 
t  «  time 

sgn®'"  slgnxim  operator  -  0  if  (O'  -  0);  ■  -1  if  5L  0) 
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5.0  CONCLUSIONS  AND  PROGRAM  FOR  NEXT  INTERVAL 


5 . 1  Practical  Implementation 

The  realization  of  equation  (13)  Is  most  easily  visualized 
as  a  real  time  analog  synthesis.  Fast  log  amplifiers  or  Schmitt 
triggers  are  required.  Our  first  attack  was  to  build  up  the  switching 
function  by  this  technique.  It  Is  now  felt  that  a  faster  than  real 
time  network  analog  of  the  linear  plant  with  a  sampling  circuit  and 
comparator  to  close  the  loop  will  be  both  economically  and  technically 
better.  This  device  will  be  decrlbed  In  Section  5.3. 

Exhaustive  catalog  searches  have  led  us  to  the  Diehl  Mfg. 
Corp.  products  as  the  most  useful  source  of  medium  power  servo  com- 


ponents . 

We  have 

in  fact  ordered: 

a. 

Diehl  No. 

FPF85-18-1 

Servo-motor 

b. 

Diehl  No. 

FPF85-29-1 

Servo-motor  w/Tach 

c. 

Diehl  No. 

FDE-15-1 

D.C.  Tachometer 

Extemal  torque  Is  to  be  provided  by  a  small  hysteresis 
synchronus  motor.  For  reasons  of  economy,  the  Inertial  load  Is  to 
be  placed  at  the  motor  shaft  with  a  precision  anti-backlash  geartraln 
(Sterling  T628,  299.7:1)  driving  the  control  transformer. 

In  order  to  maintain  a  flexible  dual  mode  boundary  It  was  deemed 
advisable  to  construct  a  full  power  saturating  amplifier.  This  Item 
Is  In  design  at  present  and  will  be  completed  by  the  motor  delivery 
date,  approximately  1  December  1961. 

5.2  Problem  Areas 

a.  Sensor  measures  Sin  e,  not  e,  this  is  demonstrably 
stable,  but  leads  to  non- optimality  for  steps  larger  than  \\/2  radians. 
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b.  Sensor  output  is  60  cycle  subcarrier,  i.e.  ,  there,  exists 
an  inherent  16  mlllisec.  time  delay  between  samples. 

c.  Component  drift  and  long  term  accurace  prohibits  use  of 
practical  multipliers  and  log  amplifiers  in  computation  of  real  time 
switching  functions. 

d.  Motor  controller  should  have  peak  power  capability  of 
600  watts  to  effect  smooth  torque  transition  across  dual  mode  boundary. 
Solid  state  switching  imposes  a  severe  restriction  at  this  power 
level . 

e.  Noise  in  tachometer,  in  input  signals  and  computer 
cannot  be  comensated  for  by  any  known  technique  other  than  construction 
pf  a  wide  dual  mode,  boundary.  These  effects  in  a  nonlinear  discontin¬ 
uous  controller  have  not  be  mathematically  investigated  and  this  Is 
beyond  the  scope  of  the  present  project. 

5.3  Proposed  Solutions 

a.  The  use  of  induction  potentiometers  would  cure  this 
difficulty,  however,  noise  in  the  difference  amplifier  would  intro¬ 
duce  as  many  problems.  This  solution  increases  the  cost  of  the  system 
without  improving  performance,  noticeably.  It  has  been  observed  that 
the  choice  of  gear  train  can  all  but  eliminate  this  ambiguity.  In 
any  event  the  problem  arises  only  for  large  values  of  e  and  causes 
only  one  ambiguous  switching  even  then. 

b.  At  runaway  velocity,  a  16  millisecond  sampling  delay 
could  cause  a  1,2°  error  in  the  output  position  of  the  device.  Clearly 
this  would  not  be  a  steady  state  error  but  would  be  a  contribution 

to  the  expected  transient.  Smoothing  and  linear  extrapolation  can  be 
accomplished  cheaply  and  should  go  a  long  way  toward  a  cure. 
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c. .  Analog  techniques  are  avoided  in  the  following  scheme: 

Construct,  a  network  analog  which  operates  faster  than 
real  time.  This  network  which  will  involve  active  devices,  is  linear 
and  simulates  the  physical  plant.  It  will  be  shown  rigorously  in 
the  next  report  that  if  the  real  time  state  variables  are  applied 
to  the  simulator  as  initial  conditions  then  an  a  priori  selection 
of  the  sgn  ct  will  enable  determination  of  the  actual  sgn ly  as  soon 
as  the  simulated  error  rate  irakes  its  .zero  crossing.  This  deter¬ 
mination  can  be  done  with  diode  logic. 

It.  appears  to  be  possible  to  achieve  good  accuracy 
at  sampling  rates  in  excess  of  400  pps.  A  detailed  description 
will  be  provided  of  the  entire  simulator- computer. 

d.  Only  rarely  will  the  full  reverse  power  of  600  watts 
be  required.  A  compromise  measure  has  been  implemented.  A 
saturating  250  watt  transistorized  controller  is  under  construction 
and  the  dual  mode  region  will  be  adapted  to  this. 

e.  The  effects  of  noise  will  not  be  studied  until  our 
computer  error  analysts  is  instituted. 

5.4  Present  Status 

a.  Hathematical  and  computational  work  Is  60%  complete. 

The.  remainder  will  be  devoted  to  computer  sy-.-ithesis. 

b.  Design  of  electronic  ."omponents  is  20%  complete. 

c.  Mechanical  design,  block  sy.ithesis  .and  component 
selection  is  90%  complete 

The.  remaining  40%  of  the  mathematical  and  computational 
work  will  be  in  the  area  of  error  analysis.  .A  computer  routine 
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has  been  obtained  which  will  enable  us  to  construct  the  topology 
of  the  dual  mode  region  in  order  to  prevent  extraneous  switching 
of  the  servomechanism  in  the  presence  of  random  additive  noise  and 
errors  in  the  calculation  of  the  true  switching  curve.  Computa¬ 
tional  work,  needs  to  be  done  in  the  calculation  of  the  linear  feed¬ 
back  function  to  be  used  within  the  dual  mode  region.  This  work 
is  essenti'  lly  st3te-o.f 'the-art  control  systems  engineering  and 
can  be  done  relatively  quickly. 

Those  equipments  which  have  already  been  designed  are  the 
error  sampling  device,  the  linear  portions  of  the  switching  compu¬ 
ter  and  part  of  the  high  power  control  amplifier.  Design  work 
remains  to  be  done  on  the  external  torque  drive  source,  the  non¬ 
linear  portions  of  the  motor  control  conputer  and  the  high  power 
stage  of  the  linear  amplifier. 

The  remaining  mechanical  design  that  is  to  be  completed  is  the 
establishment  of  error  criterion  and  tolerances  for  the  machanical 
components  associated  with  the  gear  train  and  external  torque 
drive  motor. 

Test  and  evaluation  will  be  performed,  beginning  as  soon  as  the 
device  has  been  fabricated.  It  Is  expected  that  two  man-months 
of  work  will  be  involved. 

5.5  Summa  ry 

In  summ-iry,  we  feel  that  great  progress  has  been  made 
in  establishing  s  synthesis  technique  which  will  be  useful  for 
engineers.  The  mathematical  work  has  mcovered  new  but  apparently 
soluble  problem.s  connected  with  the  mul tipli-connected  phase-space. 
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The  engineering  work  has  disclosed  a  new  but  relatively  simple 
mechanization  of  che  switching  function  using  sampled  data 
techniques.  Several  practical  problems  have  appeared  in  the 
project  and  have  been  isolated  and.  either  minimi.zed  or  solved. 

The  final  device  w’hen  delivered  will  represent  a  new  step  forward 
in  the  art  of  relay  actuated  control  systems  as  well  as  pointing 
out  several  new  mathematical  difficulties  and.  therefore  should  be 
a  significant  development  in  control  system  technology. 
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